Scaling theory of magneto-resistance in disordered local moment ferromagnets 
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We present a scaling theory of magneto-transport in Anderson-localized disordered ferromagnets. 
Within our framework a pronounced magnetic-field-sensitive resistance peak emerges naturally for 
temperatures near the magnetic phase transition. We find that the resistance anomaly is a direct 
consequence of the change in localization length caused by the magnetic transition. For increasing 
values of the external magnetic field, the resistance peak is gradually depleted and pushed towards 
higher temperatures. Our results are in good agreement with magneto-resistance measurements on 
a variety of disordered magnets. 



Complex magnetic materials showing strong magneto- 
resistance have simultaneously been at the focus of the 
attention of the magnetic recording industry and the 
field of strongly correlated electron systems. As a con- 
sequence, the interplay between electronic and magnetic 
degrees of freedom has been intensely studied both exper- 
imentally and theoretically during the last four decades. 

One of the first theoretical frameworks to provide 
guidance for understanding magneto-resistance measure- 
ments on magnetic metals were provided by de Genncs 
and Friedel |l| who predicted that long range magnetic 
fluctuations near the Curie temperature Tc will result in 
a cusp- like singularity of the resistivity p(T). As subse- 
quent experiments failed to observe such singular behav- 
ior, Fisher and Langer 0] pointed out that p(T) shows 
no anomaly, but it is dp(T)/dT that bears resemblance 
to a cusp near T c , and it is caused by short (not long) 
range magnetic fluctuations. 

While the work of Fisher and Langer provided for 
several decades a very useful theoretical framework to 
understand magneto-resistance in itinerant magnets, re- 
cently discovered complex magnetic materials, such as 
manganites and diluted magnetic semiconductors 
do show a resistivity peak directly in p(T) and a corre- 
sponding anomaly in a relatively large temperature win- 
dow near T c . For magnetic semiconductors and for some 
of the manganites the before-mentioned resistivity peak 
appears to be located precisely at T c . As pointed out 
by Littlewood and his coworkers 

13. 

such behavior, espe- 
cially for non-metallic samples, seems to fall outside the 
range of applicability of the Fisher-Langer model. Sev- 
eral proposals have been made, involving magnetic po- 
larons ||, critical spin- flip scattering 0|, thermal mag- 
netic fluctuations within the six-band model Q or Na- 
gaev's magneto-impurity scattering model Q . Most of 



these theories focus on metallic samples and on produc- 
ing the resistive anomaly near Tc, while ignoring non- 
metallic samples and at least part of the full temperature 
range. 

The available experimental results, however, provide 
severe constraints if one requires the theoretical frame- 
work to reproduce not only the relative shape and po- 
sition of the anomaly in a narrow temperature window 
near Tc (fluctuation regime), but the actual magnitude 
of the anomaly, together with the resistive behavior in 
the paramagnetic and the magnetically ordered phase as 
well. Furthermore, the theory should explain the fact 
that, with the application of an external magnetic held, 
the experimentally observed anomaly is simultaneously 
reduced in height and pushed towards higher tempera- 
tures and the magneto-resistance curves for different ex- 
ternal fields never cross. A successful theoretical frame- 
work for the magneto-resistance anomaly in disordered 
magnets should simultaneously satisfy all the experimen- 
tal constraints mentioned above. 

In this paper we show how all properties described 
above can be explained as a consequence of the in- 
terplay between disorder induced localization transition 
and magnetism, without invoking additional mechanisms 
such as Jahn- Teller distortion, probably important for 
some manganites [2^ . Our theory is based on simple 
scaling argument and is relevant for a large number of 
systems including some of the manganites and a number 
of magnetic semiconductors. 

Let us first discuss how the original one-parameter 
scaling theory of localization has to be modified to de- 
scribe magnetic systems. In the standard scaling the- 
ory of localization one argues that the typical dimension- 
less conductance g(2L) of a cube of linear size size 2L 
is uniquely determined by the typical dimensionless con- 
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ductance g(L) of its pieces of size L. This assumption is 
summarized in the following scaling equation: 

The (3(g) function above depends exclusively on g, the di- 
mension, and the symmetry of the Hamiltonian, however, 
it does not depend on other microscopic properties of the 
disorder. While (3(g) can be evaluated for large values of 
g perturbatively, in order to determine its full shape, nu- 
merical computations are needed. For very large g, (3 « 1 
in three dimensions, while for very small g it scales as 

/?(<?)« in 5 mm 

The scaling equation above can be used as a powerful 
tool to analyze the localization transition. Very impor- 
tantly, in three dimensions the (3 function always van- 
ishes at a critical value g c of g, (3(g c ) = 0. This zero of 
the (3 function is associated with the localization transi- 
tion: Consider a small piece of mesoscopic size Iq having 
a microscopic conductance go. If go > g c , then the zero 
temperature dimensionless conductance increases as we 
increase L and asymptotically behaves as g(L) ~ crL d_1 , 
with cr(<?o) the conductivity. For g < g c , on the other 
hand, the conductance scales as g(L) ~ exp(— 2L/£), 
with £(go) the localization length. From it follows 
that for so ~~ > ffc the localization length diverges on the 
localized side as £ ~ (g c — go)~ v while the conductiv- 
ity goes to on the metallic side of the transition as 
a ~ (go ~ 9c) v m three dimensions. The critical expo- 
nent v is related to the slope of the beta function at g c 
as \ jv = df3/d\ng\ 3c . 

Let us now consider a disordered local-moment fer- 
romagnet, in which local moments Hi coupled to some 
charge carriers are responsible for the magnetism. The 
first observation we make is that - compared to electronic 
processes - spin fluctuations are usually slow. This is es- 
pecially true in the vicinity of the ferromagnetic phase 
transition where the conductance peak of our interest 
appears. This implies that for transport properties the 
magnetic moments can be treated as static scatterers, 
and can be replaced by classical spins (li. While scatter- 
ing from these static magnetic moments itself is typically 
not sufficient to lead to localized charge carriers, it can 
substantially increase the effect of static disorder, and 
help to localize charge carriers. 

The T = conductance of a sample of size L will thus 
depend on the par ticular distribution of magnetic mo- 
ments P({^}) One can argue, however, that the 
13 function JQl should depend on this distribution only 
through the conductance go. Thus the effect of mag- 
netic moments appears in two ways (a) It determines the 
appropriate symmetry class of the [3 function (Unitary 
Ensemble = UE), and (b) enters the microscopic conduc- 
tance go — go(P{uj}) = go[P]- In general, go is therefore 
a complicated non-universal function of temperature, and 
magnetic field, 50 = go(T/Tc, H/Tq) which also depends 



on the microscopic details of the system. However, once 
one knows this function and the /3-function, one can use 
the scaling theory to determine transport properties of 
the system, as we shall discuss below. One possibility to 
determine go is to perform, e.g., a Monte Carlo simula- 
tion for a small system, and use the computed conduc- 
tivity as go- 

The scaling equation provides us the conductivity (or 
localization length) of the electrons, provided that the 
electronic temperature is zero. In a real system, how- 
ever, finite temperature has a dual role: On one hand, 
it changes the distribution P and thus the value of 
go{P), but one must also take it into account through 
the temperature of the conduction electrons. This is a 
rather complicated task on the metallic side, where the 
temperature-dependent dephasing length L max = L V (T) 
provides a cut-off for the scaling, with L V (T) being a 
non-universal function of the temperature T. It is, how- 
ever, simple to incorporate the effects of finite electronic 
temperature on the insulating side via Mott's variable 
range hopping formula |2(ij : 

^ex P {-C(M^) 1/4 }, (2) 

where A 5 (P) = l/N £ d (P), with N the density of states 
at the Fermi level, and C is a constant of the order of 
unity 14]. Note that in this formula £ must be deter- 
mined from the integration of the scaling equation, and 
it therefore depends on go(P) and thus on the tempera- 
ture through the distribution P. 

In order to illustrate the ideas described above, let us 
consider the disordered Kondo lattice with classical spins: 

(i,j),a i,a 
-\- J ^ ^ Qic\ a O~ aa' C-ict' ■ (3) 

Here t = 1 denotes the hopping amplitude of the con- 
duction electrons on a lattice, c\ a creates an electron 
with spin a at site i, and J is the exchange coupling 
between the spin of the electrons and the classical lo- 
cal moments flj. The e^-s in Eq. J3J denote random on- 
site energies, which we generate with a uniform distribu- 
tion between ±W. In what follows, we shall concentrate 
on the J ^> W,t limit of Eq. which is relevant for 
strongly spin-polarized systems such as manganites [2^ 
and some ferro mag netic semiconductors with polarized 
impurity bands [ljfl . In this limit Eq. © simplifies to: 

H J=00 = - ^T Ujala, + ^ eiajai , (4) 

where the a, denote spinless Fermions corresponding to 
the original Fermions aligned antiparallel with the local 
moments, and iy = e~ lipij (1 + fljfIj)/2, with ipij a Berry 



3 



phase that depends on the directions Qi and Qj |15| . 
The electronic properties of Eq. (0J have been analyzed 
for completely radom spin orientations in Ref . . Here, 
however, we want to study the effect of the ordering of 
the spins on the electronic properties. 

In principle, the finite temperature distribution func- 
tion of the spins could be computed using an effective 
action that one obtains by integrating out the conduc- 
tion electrons at temperature T. Instead of doing this, 
we shall replace this effective action by a simple mean 
field theory, and assume that the distribution of Qi is 
simply given by 

P ((U) = Jexp[fi?a] , (5) 

where a = (H + K(£l z ))/T]. Here H denotes an exter- 
nal field, K is an effective exchange coupling between the 
spins, and the magnetization m = (fl z ) must be deter- 
mined self-consistently. It is not difficult to see that in 
this case P(fli) depends in fact on a single parameter 
m = (£l z ), which is a universal function of t = T/Tc 
and h = H/Tc, with Tc — if/3 being the critical tem- 
perature. We computed m(t, h) numerically by solving 
simple transcendental equations. 

To obtain the phase diagram of the model defined by 
Eqs. J2J an d © we carried out a careful transfer matrix 
analysis within the formalism following the original work 
of McKinnon and Kramer |13l 1 1 8| . The main results of 
our analysis are summarized in Fig.^for a carrier concen- 
tration 0.5 electron/spin. Firstly, we find that larger and 
larger disorder is needed to obtain a localized phase for 
increasing m (see the phase boundary in Fig^,). In other 
words, aligning spins delocalizes electrons, and leads to a 
decrease of the resistivity. We performed a scaling analy- 
sis of the Ljapunov exponents and found that for m < 0.9 
all data collapsed to a single scaling curve, independent 
of the specific value of m, confirming the single parame- 
ter scaling hypothesis made earlier. Note, however, that 
the m ss 1 data could not be collapsed with the m < 1 
data. This is quite natural, since for m = 1 the Hamilto- 
nian belongs to a different symmetry class (Orthogonal 
Ensemble). The scaling analysis also made us possible to 
estimate the /3 function shown in Fig. \T]p. The critical 
exponent v w 1.4 is in good agreement with the earlier 
results of Ref. [lj| for the unitary ensamble. 

The phase boundary can be qualitatively understood 
if we assume that the microscopic conductance is pro- 
portional to the conductance of a single bond. In the 
most naive approximations, the effect of the magnetic 
field is just to reduce the effective value of the hopping, 
t 2 — > t 2 s = t 2 f(m), and the conductance is roughly pro- 
portional to ~ tl 



9o(m) R 



(6) 



where 



go (to = 1) denotes the microscopic con- 
ductance for fully aligned spins. The function /(to) = 



9o 
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FIG. 1: Top: Zero temperature localization phase diagram 
of the Kondo lattice model defined in the text obtained from 
a finite size scaling analysis of Ljapunov exponents (dots). 
The localization phase transition is shifted towards stronger 
disorder as we align the spins. The dashed line denotes the 
phase boundary obtained using the naive estimate, Eq. ©. 
Bottom: Beta function obtained numerically from the scaling 
analysis. The inset shows the divergence of the localization 
length and the resistivity at the critical disorder for the case 
m — 0. We estimate the critical exponent to be v ~ 1.4. 



f(m(h,t)) can be obtained in this approximation from 
the phase boundary, which is determined by the con- 
dition go(m) = g c , and is simply given by W c (m) — 
W c (m = 1)\//(to). As shown in Fig.^,, the simple func- 
tion /(to) = (1+to 2 )/2 gives a very reasonable agreement 
with the numerical phase boundary. 

Having /(to) and the /3-function at hand, we can 
now carry out the program outlined above and combine 
Eqs. Q, J2J), and JfjJ to determine the temperature- 
dependence of the conductivity in the localized phase. 
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For the sake of simplicity, we shall assume that we are 
close enough to the metal-insulator transition, and ap- 
proximate the /^-function as (3{g) ss ^\n(g/g c ) with 
v ks 1.4. Since v is not far from 1, this is a reason- 
able approximation. In this case the resistivity can be 
expressed as 



lug = A 



In 



3o (to) 



3i//4 



1 



(7) 



where A is a constant of the order of unity, and 
g c /go(Tn) = B I f(m(t, h)). The constant B here mea- 
sures simply the distance from the localized phase, B = 
g c / go(m = 1), and /(to) is the scaling function in Eq. Jfjjl. 

Typical results are summarized in Fig. [3 for B = 1.5 
and B — 4. The resistivity curves are strikingly simi- 
lar to the ones measured in various compounds in or in 
the vicinity of the localized phase, and clearly display a 
lar ge p eak at Tc and a giant negative magneto-resistance 
pi. l22j. This peak is simply a consequence of reduc- 
ing the localization length while entering the magnetic 
phase, and has nothing to do with critical fluctuations 
(which may also lead to additional contributions [2l|). 
The magnetic field dependence of the data also agrees 
qualitatively with the one seen in the experiments: The 
peak is getting flat and shifts upward with increasing 
magnetic field. One of the most important properties of 
the experimental data is that the resistivity curves cor- 
responding to different magnetization do not cross. The 
theory of Ref. e.g. does not seem to satisfy this crite- 
rion [2j , while in our theory, this is a natural consequence 
of the reduction of £. Note that in the absence of localiza- 
tion effects / disorder, the resistivity would not display 
a peak at T c 0, H^. We find a similar peak structure 
in the metallic phase, however, there the precise shape 
of the anomaly depends also on the assumption made 
for the temperature dependence of the dephasing length, 
L V (T). 

In conclusion, in the present paper we have studied the 
interplay of disorder and magnetization in disordered lo- 
cal moment ferromagnets. We proposed a unified frame- 
work to study the localization phase transition in these 
materials and argued that a unique beta function can 
be used to describe the localization phase transition in 
these materials. We verified the above hypothesis for a 
simple model of disordered local moment ferromagnets. 
The scaling approach of this paper allowed us to estimate 
the temperature and magnetic field dependence of the re- 
sistivity in the localized phase of the mean field model 
studied. The obtained resistivity curves display a peak in 
the resistivity at the critical point, due to the interplay 
of magnetism and disorder. This resistivity maximum is 
gradually suppressed and shifted towards higher temper- 
atures upon application of a magnetic field, and the com- 
puted resistivity curves do not cross. Our simple theory 
thus seems to explain all basic features of the resistivity 
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FIG. 2: Resistivity computed from Eq. J7J for the localized 
phase of the mean field model discussed in the text. We used 
A — 2 and v = 1.4 for the curves presented. The top fig- 
ure shows the resistivity for parameters not very far from the 
localization phase transition (B — 1.5), while the bottom fig- 
ure displays resistivity curves computed deep in the localized 
phase (B = 4). A peak appears at T c due to the interplay of 
magnetic ordering and localization, and is shifted to higher 
temperatures upon application of magnetic field. 



anomalies observed in many feromagnetic semiconduc- 
tors in the 'localized' phase and some of the manganites. 
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